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Abstract 

In this work we investigate small clusters of helium atoms using the hyper spherical harmonic 
basis. We consider systems with A = 2,3,4,5,6 atoms with an inter-particle potential which 
does not present a strong repulsion at short distances. We use an attractive gaussian potential 
that reproduces the values of the dimer binding energy, the atom-atom scattering length, and the 
effective range obtained with one of the widely used He-He interactions, the LM2M2 potential. In 
systems with more than two atoms we consider a repulsive three-body force that, by construction, 
reproduces the trimer binding energy of the LM2M2 potential. With this model, consisting in 
the sum of a two- and three-body potential, we have calculated the spectrum of clusters formed 
by four, five and six helium atoms. We have found that these systems present two bound states, 
one deep and one shallow close to the threshold fixed by the energy of the {A — l)-atom system. 
Universal relations between the energies of the excited state of the A-atom system and the ground 
state energy of the (^4 — l)-atom system are extracted as well as the ratio between the ground state 
of the A- atom system and the ground state energy of the trimer. 

PACS numbers: 31.15.xj, 31.15.xt, 36.90. +f, 34.10,+x 
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I. INTRODUCTION 



Systems composed by few helium atoms have been object of intense investigation from 
a theoretical and experimental point of view. The existence of the He-He molecule was 
experimentally established using diffraction experiments jl-4]. Its binding energy E 2 b has 
been estimated to be around 1 mK and its scattering length ao around 190 a.u. This makes 
the He-He molecule one of the biggest diatomic molecules. On the theoretical side, several 
He-He potentials have been proposed; in spite of different details and derivations, all of 
them share the common feature of a sharp repulsion below an inter-particle distance of 
approximately 5 a.u.. 

Another important characteristic of the He-He interactions is their effective range tq « 
13 a.u.. Accordingly, the ratio an/Vo is large enough (> 10) to place small helium clusters into 
the frame of Efimov physics O]. As shown by Efimov, when at least two of the two-body 
subsystems present an infinitely large scattering length (or zero binding energy) an infinite 
sequence of bound states (called Efimov states) appear in the three-body system; their 
binding energies scale in a geometrical way and they accumulate at zero energy. The scaling 
factor, e~ 2n / s ° ~ 1/515.03, is universal and depends only on the ratio between particle masses 



(fort 
Ref 



a 



iree identical bosons so ~ 1.00624), not on the details of the two-body interaction (see 
for a review). For a finite an/r ratio, the number of the Efimov states has been 
estimated to be N = (s /vr)ln|ao/ro| [6|]; for example, the (bosonic) three 4 He system presents 
an excited state just below the atom-dimer threshold that has been identified as an Efimov 
state. 

Triggered by this interesting fact, several investigations of the helium trimer have been 
produced, establishing that its excited state is indeed an Efimov-like state (see for example 
Refs. js-10]). In addition, analysis of the atom-dimer collision in the ultracold regime have 



also been performed llMl3] 



One of the main difficulties in solving the quantum mechanical problem in the case of 
three helium atoms results to be the treatment of the strong repulsion at short distances of 
the He-He potential. Specific algorithms have been developed so far to solve this problem. 



The Faddeev equation has been opportunely modified 14j . Moreover, the hyperspheri 



cal adiabatic (HA) expansion has been extensively used in this case 13| (for a review see 



Ref. (15J). However, due to the difficulties in treating the strong repulsion, few calcula- 



tions exist for systems with more than three helium atoms. For example, in Ref. 16j the 
diffusion Monte Carlo method has been used to describe the ground state of He molecules 
up to 10 atoms, and in Ref. 17j a Monte Carlo technique has been used to construct the 
lowest adiabatic potential in systems with 3 and up to 10 helium atoms. On the other hand, 
description of few-atoms systems using soft-core potentials are currently operated (see for 
example Ref. |18|). 

Therefore, the equivalence between hard- or soft-core-potential descriptions needs some 
creation. In a recent work Q, an attractive He-He gaussian potential has been used 
to investigate the three 4 He system. In absence of direct experimental data, the two-body 
potential has been designed to reproduce the helium dimer binding energy E2b, the He-He 
scattering length ao, and the effective range r$ of the LM2M2 potential 20]. This gaussian 
potential can be considered as a regularized-two-body contact term in an Effective Field 
Theory (EFT) approximation of the physics driven by the LM2M2 potential 2l|. It should 
be noticed that two potentials predicting similar values of ao and r predict similar phase 
shifts in the low energy limit and, therefore, even if their shape is completely different, they 
describe in an equivalent way the physical processes in that limit 2l|. The equivalence is 
lost as the energy is increased, when the details of the potential become more and more 
important. 

Extending the study to the three-body system, differences between the attractive gaussian 
and the LM2M2 potentials are immediately observed. For example the trimer ground state 
energies differ by more than 15% (see Table H]). A natural way to restore the equivalence 
between the two potentials is by the addition of a three-body soft-term force to the gaussian 
potential. On the other hand in an EFT treatment of the three-boson system with large 
scattering length, a three-body-contact term is required at leading order (LO). Its strength is 
usually determined by fixing a three-body observable as for example one of the trimer bound 
state energies. After this choice cut-off independent results are obtained 22j . Following this 



ideas, and based on Ref. 



19], in the present work we have considered a gaussian-hypercentral 



three-body force with the strength fixed to reproduce the LM2M2 ground state binding 
energy of the three-atom system. The quality of this description has been studied for 
different ranges of the three-body force. 

Using the two-atom and three-atom systems to fix the model interaction, we have analyzed 
heavier systems, up to A = 6 atoms. The numerical calculations are performed by means of 



the hyperspherical harmonic (HH) expansion method with the technique developed recently 



by the authors in Ref. [23|. In this approach, the authors use the HH basis, without a 
previous symmetrization procedure, to describe bound states in systems up to six particles. 
The method is based in a particular representation of the Hamiltonian matrix, as a sum 
of products of sparse matrices, well suited for a numerical implementation. Converged 
results for different eigenvalues, with the corresponding eigenvectors belonging to different 
symmetries, have been obtained. As a novelty, in the present work we extend the formalism 
to treat a three-body force. Moreover, as we are dealing with atoms of 4 He, only the the 
spectrum corresponding to totally symmetric eigenstates are of interest. 

After fixing the strength of the three-body force to correctly describe the LM2M2 three- 
body ground state = 126.4 mK, we have calculated the first three levels of the spectrum 
with total angular momentum L = of the A = 4, 5, 6 systems. In the three cases we have 
found that the first two levels are bosonic bound states, one deep, E<®, and one very shallow, 
E^b ) close to the threshold formed by the A — 1 system plus one atom. The third state in 
all cases belongs to a mixed symmetry with an energy above the threshold and therefore 
not representing a bound state. The appearance of only two bound states in this systems 



is in agreement to previous calculations 171 ] ■ This fact has been observed in A = 4 and 



interpreted as a consequence of the Efimov like spectrum of the A = 3 system [24( . It should 
be noticed that, whereas converged results can be found in the literature for the ground 
state of the many atom systems, the energy of the excited states is much more difficult to 
calculate and only rough estimates are available. 

To gain insight on the shallow state, we have varied the range of the three-body force 
(maintaining fixed the three-body ground state energy) and we have studied the effects of 
that variation in the A = 4,5,6 spectrum. In the range considered, the variation produces 
small effects in the eigenvalues; however, it is crucial to determine if the shallow state is 
bound or not with respect to the A — 1 threshold. Interestingly, we have observed that 
when the ranges of the two- and three-body forces have a ratio of about a/2, the ratio 
between the shallow- and ground-state energy is E^/E^_ lb ~ 1.01 — 1.02, in agreement 



with Refs. 



25 



26] . This analysis confirms previous observations that each Efimov state in 
the A = 3 system produces two bound states in the A = 4 system. Furthermore, we have 
found El$JE$ « 4.5, E$/E$ « 10.5, and E$/E$ « 18.5, which is in agreement with 



Refs. 



26 



27] 



The paper is organized as follows. In the next Section II we describe the two- and three- 
body forces we used in our calculations to reproduce LM2M2 data. In Section III the results 
for the bound states of the A = 3,4,5,6 He clusters are collected whereas the conclusions 
are given in the last section. Some technical details of the method have been summarized 
in the Appendix. 



II. SOFT-CORE TWO- AND THREE-BODY HELIUM POTENTIAL 

As mentioned in the Introduction, the 4 He- 4 He interaction presents a strong repulsion at 
short distances, below 5 a.u. This characteristic makes it difficult a detailed description of the 
system with more than four atoms. Accordingly, in the present work we have studied small 
clusters of helium interacting through a soft-core two- and three-body potentials which can 
be interpreted as regularized two- and three-body contact terms in a LO-EFT approximation 
of LM2M2. 

Following Refs. [lo|, Q we use the gaussian two-body potential 

V(r) = V e~ r2 / R2 , (1) 

with V = -1.227 K and R = 10.03 a.u.. In the following we use h 2 /m = 43.281307 (a.u.) 2 K. 
This parametrization of the two-body potential approximately reproduces the dimer binding 
energy E2&, the atom-atom scattering length ao and the effective range ro given by the 
LM2M2 potential. Specifically, the results for the gaussian potential are E^b = —1.296 mK, 
a = 189.95 a.u. and r = 13.85 a.u., to be compared to the corresponding LM2M2 values 
= —1.302 mK, ao = 189.05 a.u. and ro = 13.84 a.u.. As shown in Ref. [jj]], the use 
of the gaussian potential in the three-atom system produces a ground state binding energy 
= 150.4 mK, which is appreciable bigger than the LM2M2 helium trimer ground state 
binding energy of 126.4 mK. A smaller difference, though still appreciable, is observed in 
the first excited state (see Tabled]). 

In order to have a closer description to the A = 3 system obtained with the LM2M2 
potential, we introduce the following three-body interaction 

W(p ijk ) = Wo e" 2 4*/^ , (2) 

where = |(r^ + r 2 - k + r^) is the three-body hyperradius in terms of the distances of 
the three interacting particles. Moreover, the strength Wq is fixed to reproduce the LM2M2 



helium trimer binding energy of 126.4 mK. We have studied different cases by varying the 
range of the three-body force po between 4 and 16 a.u.. Specific cases with the corresponding 
results in the A = 3 system are shown in Table HI In the first two rows of the table we report 
the ground- and excited-binding energies of the trimer, both for LM2M2 potential and its 
gaussian representation. The excess of binding is evident for this last model. Successively, 
we report, for selected values of Wo and po, the binding energies obtained summing to the 
gaussian potential the (repulsive) three-body force. By construction the ground state has 
been fixed to the LM2M2 value and, in addition, we can observe that the excited state 
is now closer to the corresponding LM2M2 result, showing an extremely small variation with 
p ; the difference between the extremal values obtained for p = 4 a.u. and p = 16 a.u. is 
less than 1%. 



It should be noticed that the ranges R of the two-body force and po of the three-body 
force are somehow related. The gaussian two-body force can be thought as originating from 
a contact interaction regularized using a gaussian cutoff with A = (see for example 
Ref. [3]). This means that configurations, in the A = 2 system, in which two atoms have a 
relative momentum k > A remain outside the present description or, in other words, details 
of the interaction for distances below a/ (r 2 ) = R/V2 are not accessible. In the A > 2 
systems, three atoms interact through the three-body force when they happen to be inside 
a sphere of radius po/y/2 at the same time. It is clear that, as no information is introduced 
in the two-body system for distances below R/ \/2, from the relation p 2 - k = |(r?- + r 2 k + r ki ) 
and putting each distance at the value R/y/2, we obtain for the three-body range p 2 = R 2 . 
Since R has been fixed in order to describe two-body quantities, in the description of systems 
with A > 3 atoms we consider different values of po with particular attention at the region 
po « V2R. 



The calculations for the A > 3 systems, up to six atoms, are performed using the un- 
symmetrized HH basis. The method has been recently used to describe up to six nucleons 



interacting through a central potential (23| . A brief description of the method is given in the 
Appendix. The novelty presented here regards the implementation of the three-body force. 
Using a particular rotation of the HH basis it is possible to construct the potential energy 
as a product of sparse matrices. The Hamiltonian matrix is obtained using the following 
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orthonormal basis 



( P n\m[K\) = (/3 (a+1)/2 



mi 



Pp/2 \ -\-jLM 



(3) 



(a + m)\ 

where L™' (/3p) is a Laguerre polynomial with a = 3N — 1 and f3 a variational non-linear 
parameter. The matrix elements of the Hamiltonian are obtained after integrations in the 
p, Q spaces. They depend on the indices m,m' and [K], [K'\ as follows 



E 

Kj 



(m' [K']\H\m[K]) 

Er}ij,LM *yij,LM -irm,m' 
°[K] [K"]°[K'"][K'] V [K"][1 



h 2 /3 2 

m 

E 

i<j<k 



\K"\ \K" 



{T { rSm ~ K(K + 3N- 2)T%J5 [K , ][K] 



(2) 



Esp*ijk,LM sp*ijk,LM iirm,m' 
U [K\ \K"\ U [K"'\ [K'\ VV [K"] [K'' 

[K"][K">] 



(4) 

The matrices and have an analytical form and are given in Ref. 28(. The matrix 
elements V^^-n and W^^ K ,^ are obtained after integrating the matrices Vi2(p) and W(p) 
in p-space (we will call the corresponding matrices V±2 and W). Introducing the diagonal 
matrix D such that (\K'\ \ D \[K]) = S^k] ,\k']K(K + 3^ — 2), and the identity matrix / in 
i^-space, we can rewrite the Hamiltonian schematically as 

H = "^( (1) T ® / + ® D) + ^ B ™ > + E P^fWV^ , (5) 

i<j i<j<k 

in which the tensor product character of the kinetic energy is explicitly given. A method to 
diagonalize a matrix of this form is given in Ref. 



III. RESULTS FOR A = 4,5,6 HE CLUSTERS 

In this sections we present results for small clusters, up to A = 6, formed by atoms of 
4 He. Despite the differences observed at the level of the three-body system between the 
gaussian two-body force and the LM2M2 potential, it is of interest the computation of the 
spectrum produced by the gaussian two-body force only for the A = 4,5,6 systems. 

In Table [TT] we show the L = ground state E^l and the first two excited states E^l 

(2) 

and E A l for increasing values of the grand-angular momentum K using the unsymmetrized 
HH basis. The calculations have been performed up to K = 40 in A = 4, K = 24 in 
A = 5 and K = 22 in A = 6. It is a property of the HH basis that when all states having 
a fixed value of K are included in the expansion of the wave function, the symmetry of 

7 



the eigenvectors reflects the symmetries present in the Hamiltonian. Since the Hamiltonian 
is symmetric under exchange of the particles, the obtained eigenvectors have well defined 
particle permutation symmetry. In the present case the ground state E Ab and first excited 
state E^l of the Hamiltonian matrix for A = 4,5,6 are symmetric states and belong to 
the irreducible representation [A] with A = A. In all cases the second excited state E^l 
has mixed symmetry and belongs to the irreducible representation [A 1] with A = A — 1. 
In Table [TT] we also observe that the ground state binding energy, E A ^, has a very fast 
convergence in terms of K and can be determined with five digits; this value fixes the 
threshold of the continuum spectrum in the A + 1 system. True bound states in the A = 4 
systems are those having a binding energy bigger than the trimer binding energy of 150.4 mK 
and, looking at the table, bound states in the A = 5, 6 systems appears below the threshold 

(2) 

of 751.38 mK and 1945.2 mK respectively. Since in all cases the second excited state E A l 
results to be above the threshold, only two bosonic states are bound in the A = 4,5,6 
systems, one deep and one shallow close to the A — 1 threshold. The next bosonic state 

(2) 

appears above E A l and, therefore, it is not bound. This result confirms previous analysis in 
the four body sector that the lower Efimov state in the A = 3 system produces two bound 
states, one deep and one shallow. Here, we have extended this observation up to the A = 6 
system. The convergence of the E Ab is much slower than for the ground state, however with 
the extended based used it has been determined with an accuracy well below 1%. 

For this atom- atom potential the ratio r^ja ~ 1/14 and therefore we are not too far 
from the unitary limit, and we can make predictions for the universal ratios E^ Ah j 'E^ A _^ b 
and E^/E^. From the table we can observe that E$/E$ = 1.085, E$/E$ = 1.10 
and E$/E$ = 1.13. These results are not so close to the universal ratio of around 1.01 
indicating that effective range corrections are important. For the ratios with respect to the 
trimer ground state we have, E$/E$ = 5.01, E$/E$ = 12.97 and E$/E$ = 25.4. As 
we will see, these ratios are substantially modified when a three-body force is included. 

Now, we consider the model with both two- and three-body interaction. The pattern of 
convergence for the ground and excited states of the A = 4, 5, 6 helium systems, using the 
gaussian two-body potential plus the repulsive three-body potential with po = 14 a.u., is 
given in Table HTT1 The maximum grand angular momentum considered is K = 40 for A = 4, 
K = 24 for A = 5 and K = 22 for A = 6. As in the case in which only the two-body force 
has been considered, in all of the three cases only two bound states appears, one deep and 
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one shallow very close to the A — 1 threshold. The ground state presents a fast convergence 
with K and the accuracy can be estimate below 0.1 mK. The convergence for the excited 
state is slower and, for the values of K considered, its accuracy is given at the level of 3 mK. 
However, from the results it is well established that, with the value of po considered, the 
excited state, E^, is bound with respect to the A — 1 threshold. In fact, for A = 4 the 
3+1 threshold appears at 126.4 mK and the upper bound estimate for this state is 129 
mK. Its ratio E$/E$ is 1.020. For A = 5 the 4 + 1 threshold appears at 568.8 mK and 
the upper bound estimate for the excited state is 575 mK. Its ratio E^ /E^} is 1.011. For 
A = 6 the 5 + 1 threshold appears at 1326.6 mK and the upper bound estimate for the 
excited state is 1350 mK. Its ratio E^ /E^ is 1.018. The ratio between the trimer ground 
state and the ground states of the A = 4, 5, 6 systems are E^ /E^ = 4.5, Ef^/Ef^ = 10.5 
and E^JE^ = 18.5, respectively. These ratios are in good agreement with those given in 



Refs. 



25 



- 27l| . and represent a substantial improvement with respect to the case in which 
the three-body force is not included. At the ratio r^/a under consideration the use of the 
two-body soft-core potential alone reduces the Efimov character of ground and first-exited 
states, which is recovered by including the three-body force. It is interesting to compare the 
results obtained using the soft-core representation of the LM2M2 potential with the results 
of Refs. (quoted in Table IHip obtained using the original LM2M2 interaction. For 

the ground state the agreement is around 2% for A = 4, 5 and around 1% for A — 6. The 
agreement is worst for the excited state, however the results from Ref. Q are obtained 
using approximate solutions of the adiabatic hyperspherical equations. 

The overall agreement for A = 4, 5, 6 between LM2M2 and soft potential gives a further 
indication that at the LO in an EFT approach to the Efimov physics there is no need of a 
four-body force; this is only a side observation which is not at all conclusive for the lack of 
systematic study as a function of the cut-off. 

Moreover, in the four panels of Fig. [2] we analyze modifications to the spectrum of the 
systems we have considered when different values of Wo and po are used. The results for 
A = 3 can be extracted from Table [TTJ the A = 3 ground state is stable by construction, 
and small variations are observed for E^ . As shown in Fig. EK, E^ is always below the 
2+1 threshold. For A = 4, see Fig. [2b, the excited state E^ is above the 3+1 threshold, 
and therefore not bounded for values of p Q < 7 a.u. For A = 5, Fig. [2b, and A = 6, Fig. [2H, 
the corresponding excited states are above the 4+1 and 5+1 thresholds for for values of 
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Po < 12 a.u., and p < 10 a.u., respectively. For A = 5, 6 the results for the bound state 
present a bump with the smaller binding energy around po = 10 a.u.. To sum up, the most 
sensitive property of the spectrum as a function of p is, for A = 4, 5, 6, if the is above 
or below threshold. As previously discussed, a reasonable choice is po = 14 a.u., and around 
this value all the excited states are bound. 

Finally, in Table [IV] the results for the universal ratios are shown in for values of po = 
12, 14, 16 a.u.; we observe that small variation of po do not drastically change these values. 
It should be noticed that in the present analysis the unitary limit is not completely reached 
since the ratio between the two-body effective range and scattering length is r /a ~ 1/14. 
An analysis of the universal ratios as a function of a is in progress. 



IV. CONCLUSIONS 

In this work we have attached two different problems. From one side, we have studied 
the possibility of calculating bound and excited states in a bosonic (atomic) system up to 
A = 6 using the unsymmetrized HH expansion and considering soft two- and three-body 
forces. On the other hand, the model has been constructed to approximate the description 
of small helium clusters taking as a reference the results of the LM2M2 potential. These 
two problems are related since the LM2M2 presents a strong repulsion at short distances. 
Therefore, the possibility of using a soft-core representation of the original potential has 



been analyzed in detail. In Ref. [19[ bound states and low energy scattering states of the 
trimer have been analyzed using the soft-core representation of the LM2M2. The results 
obtained in that work were encouraging in the sense that they were found to be in close 
agreement to those obtained using the original potential. 

Here, we have extended the analysis to bigger systems. Therefore, the description of such 
systems with sufficient accuracy is of the main importance. To this end, we have used a 



method recently developed in Ref. [23| in which the HH basis is used without symmetriza- 
tion of the basis states. The basis is complete and, when all basis elements are included 
up to a certain maximum value of the grand angular momentum K, the eigenvectors re- 
flect the symmetries present in the Hamiltonian. In the particular case here considered, the 
eigenvectors have well defined symmetry under particle permutation and they can be orga- 
nized as belonging to the different irreducible representations of the group of permutations 
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of A objects, Sa- This simple fact has allowed to make an important statement regarding 
the number of bosonic bound states present in the systems under study. After the direct 
diagonalization of the A-body system we have analyzed the first three states for increasing 
values of K. We have considered very extended basis, up to K = 40 for the A = 4 system 
and K = 24 (K = 22) for the A = 5 (A = 6) system. This allowed to obtain converged 
results for the first eigenvalues of the spectrum. The first two were symmetric states having 
eigenvalues with energy below the continuum threshold (fixed by the lowest bound state in 
the A — 1 system) and therefore they represent true bound states. The third state was found 
to belong to a mixed symmetry and results to be above the threshold. This was the case for 
all the systems considered (A = 4,5,6) and it means that the next bosonic state has an en- 
ergy above the mixed state and therefore it is not bound. Therefore we have unambiguously 
determined that these systems present only two bound states. 



The two bosonic bound states has been studied for different values of the three-body 
potential rage po- This analysis is given in Fig. [2] where the position of the excited state 
moves from unbound to bound as po increases. The particular case po ~ \f2R is explicitly 
given in Table HT1 showing that in fact the excited state is slightly bound. Moreover, since the 
He-He potential predicts a large two-body scattering length, we have studied the universal 



ratios E^/E^andE^/E i 



(i) 



case (see Refs. 



^A-i)f These ratios have been studied in detail in the A 



24|-|26l|). Estimates have also been obtained for bigger systems 



4 



27] . Our 



calculations, obtained for one particular value of the ratio r^/a, are in agreement with those 
references on the universal character of these states in A > 3 systems. An analysis of the 
universal ratios as a — > oo is at present under way. 



Finally, we would like to discuss the quality of the description using the two- and three- 
body soft-core-potential model. We observe a substantial good agreement, at the level of 2% 
or better, for the ground states of the A-atom systems in comparison to the results of the 



LM2M2 potential given by Lewerenz (Ref. 16j). The excited states have been calculated in 



Ref. [17J though using a reduced Hilbert space. Comparing to those results we observe an 
agreement around 5%. From this analysis we can conclude that a four-body force will have 
effects beyond this level of accuracy. A deeper analysis in this subject is in progress. 
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V. APPENDIX 



Following Refs. 23|, |28|, we present a brief overview of the properties of the HH basis 
and its implementation without generating basis elements with well defined permutational 
symmetry. This approach allows to avoid the complications of symmetry-adapted-basis 



23 



29) 



construction, and to easily treat permutational-symmetry-breaking terms 

We start with the following definition of the Jacobi coordinates for an equal mass A body 
system with Cartesian coordinates ri . . . r_4 



X7V-7+I 



2j 



:r J -+i-X i ), j = l,...,N. 



J + 1 

For a given set of Jacobi coordinates Xi, . . . , Xjv, we can introduce the hyperradius p 

N v 1/2 / A v 1/2 / „ A v 1/2 



\ i=1 / \ i= i / V j>j / 



(6) 



(7) 



and the hyperangular coordinates Q 



N 



tt N = (Xi, . . . , X N , 02, . . . , (f) N ) , 



(8) 



with the hyperangles 4>% defined via 



X jv = p COS 07V 

£jv-i = psin0/v cos0/v„i 



= p sm 07v • • ■ sin 0j + i cos < 



(9) 



X2 = P Sin 07V • • ■ sm 3 COS 02 

X\ = p sin 07v • • • sin 3 sin 2 . 
The explicit expression for the HH functions, having well defined values of LM, is 



iV; 



JV 



J=2 



^ 1 (^i)<8)^ 2 (x 2 )U 2 . . ■®Y 1n _ 1 (x N - 1 )\ Ln _ 1 ®Y 1n (x n ) 



LM 

(10) 



with the indicated coupling scheme. The hyperspherical polynomial is 



V R 3 (<pj)=Mn 3 3 1 '<•<>-" l' 



i 0j )' J (sin Oj) K ' 'I'',!'"' 1 (cos 24> j) 



(11) 
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The set of quantum numbers [K] includes the n 2 . . .n^ indices of the Jacobi polynomials, 
the l\ . . . In angular momenta of the particles and the intermediate couplings L2 . . . Ljv-i- 
The Kj quantum numbers are defined as 

3 



EG 



2rn 



ni = 



K = K 



N 



(12) 



i=i 



K = Kn is known as the grand angular momentum, and W" /3 is a normalization factor. For 
the definition of the a a , where a can be either an angular momentum lj or a quantum number 
Kj, one needs to introduce the hyperspherical-binary-tree structure 30|. For example the 
tree of Fig. [T] corresponds to the choice of hyperangles given by Eq. Q, in which the 
coefficients specializes to = Kj + 3j/2 — 1 and = lj + 1/2. 

Hyperspherical functions constructed using different hyperspherical-coordinate definitions 



can be related using the T-coefficients 
following tree structures 



31 



321 ] . Schematically, these coefficients relate the 




' n i -in i _ 1 Ki 

ni_i=0 




(13) 



Here K{ = Ki-\ + U + 2rij = Ki_\ + U + 2rii The explicit definition of the coefficients is given 



in Ref. 



23] . Let us call y^(Q N ) the HH basis element constructed in terms of a set of 



Jacobi coordinates in which the i-th and z+l-th Jacobi vectors results from the transposition 
between particles j, j + 1 



1 

■-Xj + 
J 



y/(j + l) 2 - 2g 
j 



X i+ 1 



X; 



+ 1 



y/jj + l) 2 ~ 2(i + 1) 1 

: Xj -\- — Xj + x 

J 3 



(14) 



with all the other vectors equal to the original ones (transposed basis). The coefficients 



A 



i,LM 
[K][K'\ 



(15) 



are the matrix elements of a matrix A^ that allows to express the transposed HH basis 
elements in terms of the reference basis. The coefficients ^.r^ir^/i form a very sparse matrix 
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and they can be calculated analytically using the T- coupling coefficients and the Raynal- 



Revai matrix elements 



23| . A generic rotation between the reference HH basis and a basis 



in which the last Jacobi vector is defined be constructed as successive 

products of the -A^^k 1 } coefficients. Defining y^(fl 1 ^) the HH basis element constructed in 
terms of a set of Jacobi coordinates in which the TV-th Jacobi vector is defined x^y = Tj — r^, 
this coefficient can be given in the following form 

= / dn\yffi{n s )]*yffi{n%) = [A™ ■ ■ ■ A™] [K]m . (16) 

The particular values of the indices i±, . . . ,i n , labelling the matrices Af ± M , . . . , Aj 1 , depend 
on the pair The matrix 

Bi M = K"'~K"> ( 17 ) 

is written product of the sparse matrices Af s. 

We consider now the potential energy of an A-body system constructed in terms of two- 
body interactions 

V = Y ; V(i,j) . (18) 

i<j 

In terms of the HH basis it results 

E^W^WWC- ( 19 ) 

ij ij 

where the matrix elements of the matrix ^(p) are defined as 
<&(P) = (3^(0^)1^(1, 2)13^(0^)) = 

Sh,i[ ■ ■ ■ 5i N ,v N h 2 ,L> 2 ' ' ' s l n ,l' n Sk 2 ,k^ ■ ■ ■ S Kn ,k' n (20) 
x J ^^^^^^^^^'^-^^^^(pcos^^P^'^- 1 ' 



Each term of the sum in Eq. fll9l) results in a product of sparse matrices, a property which 
allows an efficient implementation of matrix- vector product. This procedure can be easily 



extended to spin-dependent potentials [33 1. 



We now consider a three body force depending on the hyperradius of a triplet of 
particles r i; rj, r^, 

= w(Pm) ■ (2i) 

i<j<k 
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The term in which k = 1,2,3 verifies pf 2 3 — x "n + x "n~i- ^ can be calculated on a 
hyperspherical-basis set relative to an non-standard hyperspherical tree with the branches 
attached to leaves and £at_i going to the same node. The transition between the two 
hyperspherical sets is given by the T-coefficients 

3$f (n*) = £ CSCr^^fe) , (22) 

njv-i 

where all the variable with the tilde refer to the non-standard tree. With this choice we 
simply have 

P123 = p cos (f) N , (23) 
and the fixed-rho matrix elements reads 



J (cos ^) 4 (sin cf> N ) 3N - 8 d<!> N V^- 1 ' ai<N - 2 (<t>N)v7 N ^ ,aKN ~ 2 (<P N )W(p cos <f, N ) 



(24) 



The three-body force matrix W(p) is extremely sparse, and it is diagonal on all quantum 
numbers but the grand- angular momentum. Finally the matrix W(p) in the standard basis 
is obtained by means of the T-coefficients 

Of^ N )\wip)\y^(n N )) = £ r^_- a ^ 

(25) 

In order to calculate the other terms of the three-body force, we use the matrices A^ AI , 
defined in Eq. (TIB]) , that transpose particles; with a suitable product of these sparse matrices 

^ / = < M ---< M , (26) 

we can permute the particles in such a way that = — r,-, and x./v-i = 2/v / 3(r/ — (r, + 
Tj)/2), and pf- k = x 2 N _ x + x 2 N , and the total three-body force reads 

V M = J2im i ] t W(p)V^. (27) 

i<j<k 
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TABLE I. The ground state E^ b and the excited state E^ b of the helium trimer calculated with 
the LM2M2 potential, with its gaussian representation and with the gaussian potential plus the 
three-body potential. In this case the two parameters, the strength Wq and the range po are given. 



potential 


4? ( mK ) E 


S } (mK) 


LM2M2 [11] 


-126.4 


-2.265 


gaussian 


-150.4 


-2.467 


(Wn (K) On fa u )) 






(306.9,4) 


-126.4 


-2.283 


(18.314,6) 


-126.4 


-2.287 


(4.0114,8) 


-126.4 


-2.289 


(1.4742,10) 


-126.4 


-2.292 


(0.721,12) 


-126.4 


-2.295 


(0.422,14) 


-126.4 


-2.299 


(0.279,16) 


-126.4 


-2.302 
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TABLE II. A = 4,5,6 binding energies of the ground, E\^, and the first two excited states, 
E^l and E^t) , for increasing values of the grand-angular-quantum number K using the two-body 
soft-core gaussian potential. We also report the symmetry of the states; the ground, E^), and 
the first-exited, E^, states are totally symmetric; the second-exited state belongs to a mixed 
representation. 



K 


E^(mK) Eg\mK) E®(mK) 


E$(mK) E$(mK) E®(mK) 


E$(mK) E$(mK) E$(mK) 




[4] 


[4] 


[3 1] 


[5] 


[5] 


[4 1] 


[6] 


[6] 


[5 1] 





725.98 


31.688 




1913.0 


642.84 




3773.1 


2010.7 




2 


725.98 


31.688 




1913.0 


642.84 


314.15 


3773.1 


2010.9 


1626.5 


4 


746.45 


77.971 




1941.2 


746.01 


400.95 


3807.6 


2140.1 


1719.3 


6 


750.15 


107.63 




1944.1 


778.79 


516.60 


3809.9 


2166.2 


1840.5 


8 


751.06 


124.48 


2.5177 


1945.0 


802.47 


571.03 


3810.8 


2188.6 


1882.5 


10 


751.28 


135.94 


29.401 


1945.2 


813.88 


608.58 


3810.9 


2196.4 


1909.0 


12 


751.35 


144.17 


50.336 


1945.2 


820.87 


634.25 


3810.9 


2200.8 


1923.4 


14 


751.37 


149.30 


66.672 


1945.2 


824.84 


653.19 


3810.9 


2202.7 


1931.9 


16 


751.37 


152.98 


79.082 


1945.2 


827.23 


657.59 


3810.9 


2203.6 


1936.5 


18 


751.38 


155.54 


89.069 


1945.2 


828.67 


678.86 


3810.9 


2204.0 


1938.7 


20 


751.38 


157.43 


97.021 


1945.2 


829.58 


687.87 


3810.9 


2204.1 


1939.7 


22 


751.38 


158.76 


103.54 


1945.2 


830.15 


695.23 


3810.9 


2204.2 


1940.0 


24 


751.38 


159.77 


108.90 


1945.2 


830.50 


701.29 








26 


751.38 


160.53 


113.38 














28 


751.38 


161.10 


117.16 














30 


751.38 


161.54 


120.37 














32 


751.38 


161.89 


123.13 














34 


751.38 


162.15 


125.52 














36 


751.38 


162.37 


127.60 














38 


751.38 


162.53 


129.42 














40 


751.38 


162.67 


131.02 
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TABLE III. A = 4,5,6 binding energies of the ground, E\J, and first-excited, E\i states for 
increasing values of the grand-angular-quantum number K. The three-body force parameters are 
po = 14 a.u. and Wq = 0.422 K. We also report the symmetry of the states; both the ground, E)?), 
and the first-exited, eQ) , states are totally symmetric. 



K 


<K) 4f(mK) 


E^imK) E 




E^(mK) B«(mK) 




[4] 


[4] 


[5] 


[5] 


[6] 


[6] 





538.93 


4.557 


1288.1 


365.1 


2293.8 


1109.9 


2 


538.93 


4.557 


1288.1 


365.1 


2293.8 


1109.9 


4 


561.69 


40.29 


1319.6 


460.4 


2331.8 


1237.3 


6 


566.68 


67.47 


1324.4 


497.6 


2336.6 


1273.0 


8 


568.21 


84.22 


1326.1 


527.0 


2338.4 


1307.7 


10 


568.58 


96.04 


1326.5 


542.7 


2338.7 


1323.1 


12 


568.73 


105.30 


1326.6 


554.0 


2338.8 


1334.4 


14 


568.77 


111.17 


1326.6 


561.0 


2338.9 


1340.9 


16 


568.78 


115.58 


1326.6 


565.9 


2338.9 


1345.3 


18 


568.79 


118.78 


1326.6 


569.3 


2338.9 


1348.2 


20 


568.79 


121.20 


1326.6 


571.8 


2338.9 


1350.2 


22 


568.79 


122.98 


1326.6 


573.6 


2338.9 


1351.6 


24 


568.79 


124.38 


1326.6 


574.9 






26 


568.79 


125.47 










28 


568.79 


126.33 










30 


568.79 


127.02 










32 


568.79 


127.57 










34 


568.79 


128.02 










36 


568.79 


128.40 










38 


568.79 


128.70 










40 


568.79 


128.96 










Ref.[16j 


558.4 




1302.2 




2319.4 




Ref.[17] 


559.7 


132.6 


1309.3 


597.1 


2329.4 


1346.7 
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TABLE IV. The ratios E%/E$j and E^/Ef', as a function of the three-body cutoff 



n ( an \ pW/pC) f(°)/f( ) pW/pC) pC 1 )/^ ) 



12 


4.47 


1.01 


10.33 


1.001 


18.12 


1.005 


14 


4.50 


1.02 


10.50 


1.011 


18.50 


1.018 


16 


4.54 


1.03 


10.70 


1.021 


19.06 


1.029 
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FIG. 2. Ground- and excited-state energies of the A = 2,3,4,5,6 systems as a function of pq. 
In panel a) we report the ground- and excited-state energy of A = 3 system together with the 
ground-state energy of A = 2; for all of the values of po we have considered, the excited ^4 = 3 
state is bounded. In panel b) we report the ground- and excited-state energy of A = 4 system 
together with the ground-state energy of A = 3; the excited A = 4 state is bounded for po > 7. 
In panel c) we report the ground- and excited-state energy of A = 5 system together with the 
ground-state energy of A = 4; the excited A = 5 state is bounded for po > 12. In panel d) we 
report the ground- and excited-state energy of A = 6 system together with the ground-state energy 
of A = 5; the excited A = 6 state is bounded for po > 10. 
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